Hacemos una tabla de los barrios de donde vienen esos usuarios y la distancia (al centro del barrio)
users_barrio_Ikea <- user_barrio_most[user_barrio_most$user %in% users_Ikea,]
agg_barrio_Ikea <- users_barrio_Ikea %>% group_by(bmax) %>% summarise(nn=length(user))
Juntamos estos datos con el mapa
mapb_Ikea <- mapb
mapb_Ikea <- merge(mapb_Ikea,agg_barrio_Ikea,by.x="GEOCODIGO",by.y="bmax")
Pintamos un choropleth con estos datos
mapp <- mapb_Ikea[!is.na(mapb_Ikea$nn),]
pal <- colorNumeric("magma",domain=mapp$nn)
m <- leaflet(mapp) %>% addTiles() %>%
addPolygons(fillColor= ~ pal(nn),
weight=1,fillOpacity = 0.8,
popup=~paste(DESBDT,"<br> nusers = ",nn,sep=""))